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Borrowing techniques from cosmology, I compute the power spectrum of quantum fluctuations in (2+l)-dimensional 
causal dynamical triangulations, a promising discrete path integral approach to quantum gravity. The results agree 
with those of canonical quantization to a high degree of precision, providing strong evidence for the equivalence of 
the two approaches and for the validity of the discrete method. 



1 Introduction 

Few quantum theories start their hves as quantum the- 
ories. Most are born in the famihar classical realm and 
quantized. There are many ways of performing this 
quantization, but they mostly fall into two broad cat- 
egories. The first, path integral quantization, is based 
on the principal of least action in classical mechanics. 
An integral is constructed over all field configurations 
between some initial and final boundary configurations 
with the action contributing a phase factor to each con- 
figuration. In the second method, canonical quantiza- 
tion, variables in the classical theory are promoted to 
operators on a Hilbert space. The structure of the the- 
ory is then defined through commutators of these vari- 
ables. 

When applying these techniques to general relativity 
many technical and conceptual difficulties arise. These 
problems are so severe that no one has yet been able to 
overcome them except in certain simplified cases. And 
the situation gets even worse; it is unclear if the results 
from the different quantization procedures, were we able 
to carry them out in full, would agree with each other 
[1]. So the question arises: if there are several theo- 
ries that yield the same classical behavior but describe 
very different quantum mechanics, which is the correct 
quantum theory? 

In this paper I will examine the results from two 
different approaches to quantizing gravity that can be 
carried out in full. The first is the reduced phase space 
quantization of (2+l)-dimensional gravity with spher- 
ical topology. The reduced dimensionality and simple 
topology yield a trivial classical spacetime. This simplic- 
ity is preserved when the theory is quantized, resulting 
in a quantum theory with no degrees of freedom. The 
second approach is causal dynamical triangulations, a 
lattice regularization of the gravitational path integral. 
Smooth spacetime is approximated by a triangulated 
mesh and the path integral is approximated by a sum 
over inequivalent triangulations. The sum is computed 
using Monte Carlo techniques, and the resulting space- 
time is analyzed using techniques borrowed from the 
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analysis of cosmic microwave background temperature 
fluctuations. I find that the resulting measurements re- 
veal an extremely spherical spacetime, indicating that 
the two quantization techniques yield very similar re- 
sults. 

2 (2+l)-dimensional gravity and 
canonical quantization 

The action for general relativity is the Einstein-Hilbert 
action 

Seh = y^ j^d^x^g{R-2K), (1) 

and the Euler-Lagrange equations give us the Einstein 
field equations 

R^. - \9^.R + A^M- = -^^GT^^. (2) 

As one would expect, working in 2+1 dimensions 
greatly simplifies the equations of motion. In d dimen- 
sions ([2| allows d{d — 3) local phase space degrees of 
freedom |21 13| . In d = A this gives us the 4 degrees of 
freedom that result in gravitational waves. In (i = 3, 
however, this results in zero local degrees of freedom. 
Additionally, in 2+1 dimensions the curvature tensor 
R^jLupcr depends linearly on the Ricci tensor R^y . This 
means that vacuum solutions to Q are flat if A = 
and of constant curvature if A ^ [2 , 3 . 

So it seems that choosing to work in 2+1 dimensions 
has given us a fairly boring universe to work with. It 
turns out that exactly how boring depends on the topol- 
ogy of spacetime. In general our manifolds M will de- 
compose into one-dimensional timelike intervals / and 
two-dimensional spacelike manifolds S ; more compactly 
/ X S. If we choose S to be any manifold contain- 
ing non-contractible curves (e.g., a torus) the resulting 
holonomies will yield a finite number of global degrees of 
freedom [3, 4 . For our purposes, we want the simplest 
possible scenario to work with, so we will work with a 
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spacelike manifold in which all curves are contractible. 
If we choose our constant curvature to be positive, this 
leaves us with a two-sphere S'^ . The simplest possible 
2+1 dimensional topology we can look at is then / x /S^ . 

Quantizing even this simple spacetime is far from 
straightforward. A detailed description of how one goes 
about doing this can be found in [3 . The result is that 
if one chooses to apply the classical constraints first and 
then quantize the system, one ends up with a quantized 
spacetime with no local or global degrees of freedom. 
The Hilbert space is one-dimensional; the 2-sphere that 
we started with remains a 2-sphere, around which there 
are no quantum fluctuations. 

3 Lattice quantum gravity and 
covariant quantization 

The basic idea behind quantizing a classical theory using 
the path integral approach is to form an integral over all 
possible field configurations between some fixed initial 
and final configurations. The field configurations we will 
integrate over are essentially all of the unique geome- 
tries characterized by spacetime metrics interpolating 
between some starting spatial geometry and some end- 
ing spatial geometry. Along with the Einstein-Hilbert 
action ([T]), we now have all the ingredients we need to 
write down a formal path integral for quantum gravity: 



iSEh 



(3) 



In practice, expressions like (|3| are usually calcu- 
lated using perturbation theory. This approach fails 
with gravity, however, because the perturbation series 
is non-renormalizable [5]. In order to proceed, we need 
to introduce some kind of regularization scheme. We 
will do this by approximating the smooth manifolds M 
by a collection of flat simplices glued together at their 
edges, very much like a smooth sphere is approximated 
by a geodesic dome. 

In two dimensions, the simplex used is a triangle. 
Since the triangle is flat, any curvature in the surface 
created by glued together triangles is concentrated at 
vertices v where i triangles come together, and is char- 
acterized by the deficit angle: 



8y =271 - ^^i 



(4) 



In n dimensions, triangles are replaced with n- 
dimensional simplices connected to other simplices by 
(n — 2) -dimensional vertices or "bones." The curvature 
at each bone is still described by Q , where the dihedral 
angle 6 is now measured around an (n — 2) -dimensional 



bone. The sum over deficit angles, weighted by the "vol- 
ume" Vh of each bone, is proportional to the integral of 
the scalar curvature over the surface: 



h 



(5) 



We are now ready to write down a discrete version 
of ([T]), called the Regge action [HUTj: 
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(6) 

The path integral ([3| over all unique field configura- 
tions is then taken to be a sum over all unique triangu- 
lations T weighted by Sr : 
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C{T) 



(7) 



where 



C(T) 



is a symmetry factor with C{T) being the 



order of the automorphism group for triangulation T 

3.1 Causal dynamical triangulations 

Historically, discrete approaches to gravity were con- 
structed using Euclidian building blocks. This was 
mostly done for technical reasons; the weight factors in 
^ become real, which is usually necessary for the sums 
to converge. However, these approaches were unable 
to reproduce classical spacetimes [8]. Causal dynami- 
cal triangulations (CDT) [8^ attempts to address this 
problem by taking the Lorentzian structure of space- 
time seriously from the outset. The building blocks are 
Lorentzian simplices with some (n-l)-dimensional faces 
being spacelike and some timelike. The action ([6| is con- 
structed using these building blocks and then rotated 
into the Euclidian sector to aid in computation. The 
resulting construction exhibits three distinct phases. 
Two of these phases resemble the pathological results 
from Euclidian dynamical triangulations, a crumpled 
phase and a branched polymer phase, neither of which 
seem able to reproduce classical spacetimes. The third 
phase, however, is characterized by an extended geome- 
try with small fluctuations, as one would hope for from 
a quantum theory of gravity. This phase is the result of 
building the causal structure into the theory from the 
outset, and is the main success of the CDT program. 
This paper will only focus on the (2+l)-dimensional 
construction, although the (3+l)-dimensional approach 
is very similar. Details of both cases can be found in 
[8. 

In (2+l)-dimensional CDT we use three types of 
three-dimensional blocks to build our spacetime. Each 
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block is constructed from points that lie in the constant 
time surface t = ta and the next constant time surface 
t = ta -\- I. We will label our blocks by "(number of 
points at t^, number of points at ta + 1)." With this 
convention our three blocks are the (3,1) simplex, the 
(1,3) simplex (a reflection of the (3,1) simplex) and the 
(2,2) simplex (Figure [T]). 




Figure 1: The three simplices of 2+1 CDT. 



We define the length of edges that lie in constant 



and the length of 



time surfaces to be /gpace 
edges that connect two neighboring timelike surfaces as 
^time ~ —(^a'^ ; a is an asymmetry parameter between 
timelike lengths and spacelike lengths. With these defi- 
nitions the Regge action (|6| in three dimensions becomes 



8nG 



spacelike 
^ bones b 



timehke 
bones b 



A 

SttG 



(8) 



ysimpHces s 



Following the work in [8 , Lorentzian volumes Vi and 
dihedrial angles Oi for geometrical elements i are ap- 
plied to (|8|. Then a series of relationships between all 
of the bulk variables in the manifold are used, along with 
the simplifying assumption that a = l[^ and a choice of 
units a = 1. Lastly the entire action is Wick rotated 
into the Euclidean sector. The result is [51 [TOl: 
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with 



^Even without this assumption, trigonometric identities reduce 
the action to the form|9] although the relations |10| become more 
complicated. 
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where Nq is the number of vertices in the manifold, A^3 
is the number of three-simplices, G is Newton's con- 
stant, and A is the cosmological constant. Finally, the 
discrete path integral ^ becomes: 
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(11) 



3.2 Monte Carlo Moves 



In order to numerically implement the partition func- 
tion (11), one needs to create a procedure that will ex- 



plore different allowed triangulations of a given three- 
dimensional spacetime. This is accomplished by us- 
ing the Metropolis Monte Carlo algorithm [11 with a 
series of so-called "moves," each move being a simple 
re-triangulation of a set of initial simplices. The con- 
straints on these moves are that they should: (i) pre- 
serve the existing time-slicing of the spacetime, (ii) pre- 
serve the topology of the spacelike slices and, (iii) be 
ergodic, that is, any one allowed triangulation can be 
transformed to any other allowed triangulation by re- 
peated application of these moves. In three dimensions 
there are five of these moves [8], four of which are re- 
lated by inversion symmetry. They are shown in Figure 
[2] Each move is attempted a roughly equal number of 
times and is accepted or rejected with probabilities de- 
pending on how it changes the action (|9| and how it 
effects the local geometry. The unit of simulation time 
is called a sweep and is defined as TVs attempted moves 
(i.e one attempted move for every simplex in the space- 
time). The specific CDT implementation that was used 
for this work is described in [T2l. 



4 Process 

The goal of this work is to measure the sphericity of 
our CDT quantized spacetime, specifically the ensemble 
average of the sphericity of individual spacelike slices 
across many sweeps of the simulation. There are four 
steps required to do this: (i) choose a reference spacelike 
slice in each sweep to be analyzed, (ii) embed this two- 
dimensional manifold in three dimensions, (iii) measure 
the distance from the three-dimensional center of ge- 
ometry to the embedded manifold and, (iv) analyze the 
distances across many sweeps using spherical harmonics. 

4.1 Step 1: Choose reference slices 

Although this is the least technically challenging step 
of the process, it does present some conceptual difficul- 
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in the slice) in each sweep as a reference. This ensures 
that there is at least one common bulk characteristic 
between any two given slices. See [9] for a further dis- 
cussion related to this choice and the issues surrounding 
it. 
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Figure 2: The three Monte Carlo moves that, along with their inverses, 
make up the five move set (the (4,4) move is it's own inverse). 



ties. Our (2+l)-dimensional spacetime has the topology 
S'^ xl . The timelike interval / runs from t = 1 io t = T 
for T time-slices. For computational simplicity the first 
time slice is identified with the last, changing the topol- 
ogy to S'^xS^ . When the CDT program is implemented 
as described in Section [Sj the results are a bulge of ex- 
tended spacetime and stalks with the minimal number 
of simplices required to maintain the S'^ input topol- 
ogy (Figure [3|. The extended bulge performs a random 
walk along the time axis as we execute sweeps. Since 
the time axis is a circle, there is no unique choice of 
time coordinates for the bulge, so it makes little sense 
to compare a fixed time-slice across sweeps. The solu- 
tion to this dilemma is to use the maximum volume slice 
(measured by counting the number of spacelike triangles 




Figure 3: Plots of the volume per time-slice for simulation 3. The top 
is sweep 500, before the simulation has thermalized. The bottom is 
sweep 76100 where the volume bulge is clearly visible. 



4.2 Step 2: Embed the slices in 3D 

The property of spacetimes generated by the simulation 
that is of interest is the the intrinsic curvature. This is 
defined by the deficit angle Q, which is easily measured 
for individual spacetimes. But because the number and 
location of vertices will change across sweeps, it is un- 
clear how to measure an expectation value using this 
quantity. To overcome this problem it is useful to em- 
bed the individual spacetimes in three dimensions. The 
resulting coordinate system can then be used to compare 
measurements across sweeps. Of course once the space- 
time is embedded, it is extrinsic not intrinsic curvature 
that is being measured. However, Gauss's Theorema 
Egregium states that the Gaussian curvature of a two- 
dimensional surface embedded in a flat three-manifold is 
intrinsic to the surface, and is completely determined by 
the intrinsic curvature fT3^. So in 2+1 dimensions, mea- 
suring the extrinsic curvature is equivalent to measuring 
the intrinsic curvature. 

4.2.1 The embedding algorithm 

The points that define vertices in CDT are by ne- 
cessity coordinate-free. The only information that the 
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simulation generates is the linking between them. In or- 
der to embed a two-dimensional spacelike slice in three 
dimensions for analysis, we have to choose a coordinate 
for each point based on its relationship to its neigh- 
bors. This is accomplished by feeding the logical links 
between the vertices into a force-directed graphing al- 
gorithm, specifically the spring electrical embedding al- 
gorithm [14 as implemented by Mathematica 6. 

In this graphing algorithm, neighboring vertices are 
attracted to each other by a force that is proportional 
to the Euclidean distance between them, i.e., a spring 
force. Additionally, every vertex repels every other ver- 
tex with a force that falls off with the inverse square of 
the distance, i.e., an electrical force. The total energy 
of the system is given by: 



4.4 Step 4: Spherical harmonic analysis 



X q X j 



^ 4 



neighbors 
of i 



(12) 



where Cg and Cg are constants that control the relative 
strength of the electrical and spring forces, respectively, 
and dij is the Euclidean distance between vertex i and 
j . This energy function is then minimized by iteratively 
moving each vertex in the direction of the spring force. 



The configuration of points that minimizes ( 12 ) is used 
as the final embedding of our spacelike slice. 

4.3 Step 3: Measure distances 

Given an embedding of our surface in flat three-dimensional 
space, we can, in principle, determine the extrinsic cur- 
vature from the positions r{t) of the points. In partic- 
ular, for a spherical surface, r(t) will be constant, and 
deviations from a constant value will give a measure of 
the variation of the curvature. Because of this, a mea- 
surement of r{t) for each of the spacetimes generated by 
the simulation is therefore a measurement of curvature. 
Because of this I will focus on the measurement of r{t) 
for the remainder of this paper. 

To perform this measurement, rays emanating from 
the center of geometry of a particular slice are traced 
and the point at which they intersect with the mani- 
fold is calculated. The rays used are not arbitrary, but 
rather picked so they will point to specific locations on 
a 2-sphere which correspond to different pixels in the 
HEALPix discretization scheme (more on HEALPix in 
Sec. 4.4). So what is actually measured is the distance 



from the center of geometry to the two-dimensional sur- 
face as a function of HEALPix pixel number or r{npix) . 
In order to avoid correlations that are an artifact of the 
embedding algorithm, a random rotation is performed 
about the center of geometry for each set of distance 
measurements. 




Figure 4: Examples of HEALPix pixelations at different values of 
Nside (image taken from |15p. The value Ngide is the number of 
divisions of the base pixel along one side. Starting at the top left the 
first sphere is pixelated with the base resolution pixels. The top right 
sphere is at Nside = 2. The bottom right sphere is at Nside = 4. 
The bottom left sphere is at Nside — 8- 



The measurements in Section 14.31 are taken with a 
specific method of spherical harmonic analysis in mind. 
HEALPix (Hierarchical Equal Area iso-Latitude Pixeli- 
sation) is an efficient method of pixelating a 2-sphere 
that lends itself to spherical harmonic transformations. 
It is used widely in cosmic microwave background anal- 
ysis, including data from the WMAP and Planck exper- 
iment. Details of the HEALPix implementation can be 
found in [151116]. The basic idea is that any function / 
on a sphere can be decomposed into a linear combina- 
tion of spherical harmonics Yim '• 

Imax 

/(^,0)= ^^a,„F,™(e,0). (13) 

1=0 m 

The aim can then be found by using the orthogonal- 



ity of the to invert (13): 



Jn 



(14) 



In the case of a discrete function /, evaluated at 



pixel locations {Op^<pp)^ with Npi^ pixels, (14) can be 
approximated by a sum. The approximation used by 
HEALPix is yjj 
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(15) 



p=0 
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The variance of the air, 
trum, Ci : 



Ci 



=—y 



is the angular power spec- 



(16) 



The resolution of HEALPix maps is characterized 
by the quantity Nside which represents the number of 
divisions of a base-resolution pixel [15 . Examples of 
pixelations with different values of Ngi^e are shown in 
Figure [4j 

4.5 Errors in the embedding 

Although the embedding algorithm does a very good 



job, it is not perfect. The requirement of Section [3yT1 that 
all link lengths be equal is very difficult for the algorithm 
to achieve; most links are the same length, but there are 
always some that are either bigger or smaller than the 
average. When the link lengths are different sizes, the 
geometry of our spacelike slices becomes distorted. This 
means that there is some difference between the actual 
geometry of our manifold and the geometry we measure 
in step 3 of our process. To understand how this will 
effect our measurements, we look at how each triangle 
that is intersected by a ray is distorted and apply an 
extra amount Ar to our measurements. We estimate 
Ar as follows: the change in the area of a triangle with 
changes in link length / is roughly AA ~ lAl. The 
change in the area of a spherical spacetime with changes 
in the radius r is roughly A ~ r Ar . Since the spherical 
spacetimes are constructed of triangles and the radius 
r roughly corresponds to our ray-trace measurements, 
we can equate these two quantities to get an estimate of 
Ar: Ar ~ -A/, where / is the average link length for 
a given spacetime, Al is is the difference of the average 
of a specific triangle's link lengths and the average link 
length for the spacetime, and r is the measured distance 
to the triangle in question. 

This quanity Ar is calculated for each pixel and ap- 
plied to the original measurement. Then the new dis- 
tances are are analyzed using the methods described in 
Section |4.4l The results of this additional analysis are 



included in Figures |8] and [9) which are discussed in Sec- 
tion |5l 



5 Results 

In order to test the sphericity of the output from CDT, 
five simulations were run. The values of the bare cou- 
pling constants for each of these simulations are listed 
in Table [l] A "snapshot" of the spacetime was exported 
to a ffie every 100 sweeps. The largest volume spacelike 
slice was then chosen from each of these spacetimes and 
used for analysis. Examples of several of these slices, 
embedded in three-dimensions, are shown in Figure [5] 



All of the HEALPix analysis was performed at 
^side = 32. This resolution effectively limits the power 
spectrum analysis to the ffist 100 / 's. At higher values 
than this, we start sampling the map at scales smaller 
than the map resolution, which results in noise in the 
spectrum. The power spectrum, up to / = 100, and 
HEALPix map for each of the slices in Figure [5j are 
shown in Figure |6] and Figure [7| respectively. 






Figure 5: Various max volume spacelike slices from simulation 1 em- 
bedded using the algorithm described in section [4.2| Top left is sweep 
55500, top right is sweep 71800, bottom left is sweep 92700 and bot- 
tom right sweep 99900. 




Figure 6: The HEALPix maps for the slices in Figure [5] Top left IS 
sweep 55500, top right is sweep 71800, bottom left is sweep 92700 and 
bottom right sweep 99900. 



5.1 The shape of spacetime 

As can be seen in Figure [7| individual slices vary from 
being spherical by non-trivial amounts. However, in- 
dividual slices come from particular paths in the path 
integral (|3|, and physically, only the ensemble averages 
are meaningful. It is important to choose carefully when 
picking the quantities from which to calculate these av- 
erages. Because all notion of absolute location is washed 
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Table 1: Various parameters and results for each simulation. Vinn is the total initializatio n vol ume (i.e. the number of initial simplices), T is 
the number of time-slices, ko and ks are the bare coupling constants discussed in Section |3.1| The numbers in parentheses next to the total 
samples for simulation 3 are the number of random subsamples of the full set of sweeps t hat were made. These results were then averaged to 
avoid sample bias. 63"^ is the deviation from being perfectly spherical discussed in Sectionls.ll 




Figure 7: Ci/Cq vs. I for each of the individual spacetimes in Figure 
[5] Top left is sweep 55500, top right is sweep 71800, bottom left is 
sweep 92700 and bottom right sweep 99900. 



out from sweep to sweep, averaging the distance mea- 
surements is an unreliable measure of the shape of the 
ensemble system. Each pixel will explore roughly the 
same range of allowed values, resulting in the same av- 
erage distance everywhere - the measurements will "sum 
to zero." If there were a non-spherical shape preferred 
by the system, it would be washed out in this average. 
We need a quantity that is rotationally invariant, i.e., 
some measurement that will give us the same answer for 
a given shape regardless of how that shape is oriented 
relative to a coordinate system. Luckily, the power spec- 
trum is just such a quantity. Averaging each / of the 
power spectrum across the entire ensemble will give us 
a measure of the physical shape of our spacetime. 

To test this assumption, simulation 3 was run for 
four times as many sweeps as the other simulations. 
This pool of sweeps was then sampled at different sizes. 
As sample size increases, measurements that are subject 
to the "sum to zero" effect will converge to average val- 
ues everywhere. Measurements that are free from this 



defect will remain unaffected. If the power spectrum en- 
semble average is one of these quantities, then it should 
be almost the same at any reasonably large sample size. 



Figure |lla| confirms this. The averaged power spectra 
are almost identical for every / at all sample sizes. 

Now that we have some confidence in our averaging 
procedure, we can examine the results for all of the sim- 
ulations in Table [T] The results are shown in Figure [Sj 
All of the simulations performed here seem to be very 
spherical. In each case, the / = 1 mode is only a few 
percent of the the spherical or / = mode, and higher 
modes fall off rapidly. Also, the results are very similar 
for all chosen values of ko and ks , implying that the 
shape of the spacetime in the extended phase does not 
depend on the values of the coupling constants. Finally, 
in order to give each simulation a single score to rank 
its sphericity, the following quantity is calculated: 



100 



-2 _ 1=1 



5S' = 



Co 



(17) 



SS^ is a measure of how much power is not in the 
spherical / = mode. The results for each spacetime 
are shown in Tabled] 



5.2 The fluctuations 

The spacetimes produced by CDT seem quite spherical, 
matching the results from canonical quantization. The 
other prediction of the canonical procedure is that the 
spacetime should have no fluctuations. To test whether 
the CDT spacetimes exhibit this property I will look at 
the variance of the distance measurements: 



^pix = {riripicc) ) - (r(npi^)) . 



(18) 



As with the measurement in Section [STT] it is important 
to note that (18) is not affected by the "sum to zero" 



effect, as can be seen in Figure [TTbj 

Equation (18) is calculated for every pixel in the 
HEALPix map, creating a map of the fluctuations of the 
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spacetime (Figure [TQ|). A power spectrum is then calcu- 
lated for this map. The result for each simulation can 
be seen in Figure [9j Clearly non-spherical (non / = 0) 
modes of fluctuations are very small, a result that fits 
with the predictions of the canonical quantization pro- 
cedure. There is however, power in the the / = mode 
that does not match the canonical predictions. This zero 
mode power can be attributed to lattice effects in the 
CDT model. Because the shces are discrete, when the 



largest volume time-slice is selected (see Section 4.1) it 
will typically be slightly smaller than the maximum vol- 
ume slice in the limit where the number of time-slices 
gets very large. This means that the volume I sample 
will fluctuate around its "maximum" value. Because the 
maximum volume is fluctuating and our ensemble space- 
times are very spherical, the radius of the spacetime will 
fluctuate uniformly, resulting in a large zero mode in 
the fluctuation power spectrum. To calculate an order 
of magnitude estimate for this fluctuation, the points 
marking the volumes of the time-slices around the max- 
imum volume time slice were fit to a curve using a cubic 
curve-fitting algorithm. The maximum of this curve was 
then interpreted to be the "actual" maximum volume in 
the limit where the number of time-slices becomes large. 
The average difference between the "actual" value and 
the selected value was calculated. In order to compare 
this average Ar to the zero mode of the power spec- 
trum, one must relate the Ci to the total power in a 
function / via Parseval's theorem: 



oo I 
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Combining this with (16) gives: 



„ oo 



(20) 



The function we are considering is the variance (18) 



which we will assume is constant over the sphere, and 
because it is so much bigger, the only Ci that we need 
to look at is Co , so: 



cr^ ^ Cq. 



(21) 



Additionally, since we are assuming that the space- 
time is mostly spherical r ~ W and Ar ~ V^max — 
V^measured (wc are workiug in two dimensions, so the 
volume in this case is actually an area). So we have 



Ar ^ \/Vmax - \/Vr 



measured 



1/4 



(22) 



Results of this calculation are shown in Figure [T2j 
Without any modification, the values for Ar that we 



measured are smaller than the Co s but they show simi- 
lar behavior. The fit can be made better by multiplying 
and adding a constant to Ar. These constants arise 
from the fact that the assumptions we made in deriv- 
ing (22) - constant variance and perfect sphericity - 



were not completely correct using our data set (there 
are small but non zero modes > in our simulations). 



The results in Figure 12 are an indication that the zero 
mode fluctuations are indeed induced by the CDT lat- 
tice. 
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Figure 12: The zero mode fluctuation (Cq'^^) compared to the differ- 
ence in the sampled radius from the actual radius due to lattice effects 
(Ar). 



5.3 Small Q for / > 

The power in / 's greater than zero is small in all cases, 
but it is still non-zero. If the spacetimes were truly 
spherical, with no fluctuations, all of these values would 
be zero. Is there a small but interesting dynamic in 
the theory itself that causes these modes to be active? 
Within the context of this work it is unclear where these 
modes come from. It is plausible that they arise from 
lattice effects. If so, then as the lattice spacing de- 
creases, theses modes would also decrease. Decreasing 
lattice spacing requires an increase in the initialization 
volume, though, and generating enough sweeps to calcu- 
late good statistics is problematic at high volume using 
available hardware. Studying these small, high / modes 
is therefore left to future work. 



6 Conclusions 

As Figures [S] and [9] and the analysis in Section [5] sug- 
gest, the geometry that the CDT simulations produce 
are very close to being perfectly spherical with no fluctu- 
ations. Because reduced phase space quantization leads 



9 



to a quantum spacetime with no degrees of freedom, 
these results imply that, at least in (2+1) dimensions, 
the path integral quantization implemented by CDT 
is, at least very nearly, equivalent to the canonical ap- 
proach. It remains to be seen if this result can be ex- 
tended to other 2+1 dimensional topologies. In CDT 
we have, for the first time, a tool to numerically explore 
these fascinating theories and perhaps cast some light 
on the (3+l)-dimensional case. 
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Figure 8: The relative power of modes greater than I = when tes ting the ensemble shape of the spacetime for each of the simulations 
for measurements with and without the distortion described in Section |4.5| applied. As can be seen from the plots the distortion induced by 
embedding errors is small. 
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Figure 9: The total power of ensemble spacetime shape fluctuations for each simulation with and without the distortion described in Section |4.5| 
applied. As can be seen from the plots the distortion induced by embedding errors is small. 
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Figure 10: Maps of the ensemble spacetime shape fluctuations for each simulation. 
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